# simulation setup
n_unit <- nrow(W_use) # 2430
W_sum <- apply(W_use, 1, sum)
prob_g <- mean(apply(W_use, 1, mean))

# Sensitivity parameters: overlap 
overlap_use <- c(0.0,  0.2,  0.4,  0.6)

# parameters of the outcome DGP
true_G <- 1
model_sim <- c("linear", "interactive")

# linear model
intercept_l <- 5
coef_T_l <- 2
coef_U_l <- 1.5

# interactive model
intercept <-  5
coef_T1 <- 2
coef_T2 <- 2 
coef_U1  <- 0.5
coef_U2  <- 2

# random error
sd_norm  <- 0.5